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Abstract 

In this paper, we propose a new beamforming design to maximize energy efficiency (EE) for 
multiple input single output interfering broadcast channels (lEBC). Under this model, the EE problem 
is non-convex in general due to the coupled interference and its fractional form, and thus it is difficult 
to solve the problem. Conventional algorithms which address this problem have adopted an iterative 
method for each channel realization, which requires high computational complexity. In order to reduce 
the computational complexity, we parameterize the beamforming vector by scalar parameters related 
to beam direction and power. Then, by employing asymptotic results of random matrix theory with 
this parametrization, we identify the optimal parameters to maximize the EE in the large system limit 
assuming that the number of transmit antennas and users are large with a fixed ratio. In the asymptotic 
regime, our solutions depend only on the second order channel statistics, which yields significantly 
reduced computational complexity and system overhead compared to the conventional approaches. 
Hence, the beamforming vector to maximize the EE performance can be determined with local channel 
state information and the optimized parameters. Based on the asymptotic results, the proposed scheme 
can provide insights on the average EE performance, and a simple yet efficient beamforming strategy 
is introduced for the finite system case. Numerical results confirm that the proposed scheme shows a 
negligible performance loss compared to the best result achieved by the conventional approaches even 
with small system dimensions, with much reduced system complexity. 
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I. Introduction 

A design of traditional wireless networks foeusing on high spectral efficiency has caused 
rapidly increasing energy consumption and negative impact on the environment. Therefore, 
pursuing high energy efficiency (EE) becomes an important and urgent task for future wireless 
system designs [[II. In general, the EE is defined as the ratio of the sum-rate to the total power 
consumption measured in bit/Joule. Meanwhile, coordinated beamforming schemes, which allow 
base stations (BSs) to jointly optimize their transmissions by sharing channel state information 
(CSI), are considered as a key technology in cellular networks due to its significant spectral 
efficiency improvement Q. When the EE is taken into account for a design of future wireless 
systems, cooperative transmission techniques need to be investigated with new perspectives. 

Erom the EE point of view, several papers have studied methods to maximize the performance 
various systems [|3l-[|9||. The EE maximization problem in general belongs to a class of fractional 
programming due to its fractional form, and thus is nonlinear. Nevertheless, for the special case 
of no interference among users, the problem can be transformed into an equivalent convex 
problem without loss of optimality by exploiting the pseudo concavity of the objective function 
O. As a result, the global optimal solution can be found efficiently by convex optimization 
tools GQl. Also, for the multi-user case, the same framework can be adopted by employing 
zero-forcing beamforming ffTTI . although the resultant performance is suboptimal in terms of 
EE. However, in more general scenarios with inter-user interference, the optimization problem 
for EE is non-convex, and thus it is difficult and more challenging to optimize the EE in the 
presence of interference. 

In this paper, we focus on designing a new energy efficient scheme for multi-cell multi¬ 
user downlink systems where each BS equipped with multiple antennas communicates with its 
corresponding single-antenna users. These systems can be mathematically modeled as multiple 
input single output (MISO) interfering broadcast channels (lEBC). After transforming from 
fractional programming to linear programming in Q and applying the weighted minimum 
mean square error (WMMSE) approach in ifT^ . a local optimal solution was achieved in 
lEl. However, this method requires either centralized channel knowledge or the exchange of 
additional parameters. In addition, the optimal beamforming vectors should be computed in an 
iterative manner for every channel realizations. Moreover, it is difficult to get insights on average 
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performance without resorting to Monte Carlo simulations. 

To overcome these issues, we propose a low complexity energy efficient scheme with a 
negligible performance loss compared to the best results achieved by the conventional approach. 
To this end, we first parameterize the beam vectors by the parameters associated with beam 
direction and power. With this parameterization, we then employ the asymptotic results of random 
matrix theory [[T3]l - [fT^ . More specifically, in the large system limit where the number of transmit 
antennas and users in each cell go to infinity with a fixed ratio, we identify the parameters to 
optimize the EE. Note that the beamforming vector which maximizes the EE performance is still 
constructed based on local instantaneous CSI in a finite dimension. Meanwhile, the parameters 
can be optimized by adopting the large system analysis. It is worth noting that in the asymptotic 
regime, the parameters become deterministic and the randomness according to instantaneous 
channel realizations disappears. Therefore, only second order channel statistics is required for 
the large system approach. 

In ifTTll . the sum rate (SR) maximization is performed for MIMO interference channel by 
utilizing the fact that a zero gradient value for sum rate maximization under fixed full power is 
efficiently obtained by utilizing a relationship between the SR and virtual signal-to-interference- 
plus-noise. However, this method cannot be directly applied for the EE metric since beamforming 
direction and power allocation are jointly considered for the EE maximization. Different from 
conventional EE algorithms which should be updated in each channel realization, the proposed 
method does not recalculate the parameters as long as statistical channel information remains 
constant. Thus, our long-term strategy significantly reduces complexity and the system overhead 
compared to the conventional methods. In addition, the dimensionality of the optimization 
problem is greatly reduced by virtue of the asymptotic approach. Moreover, an asymptotic 
expression of the achievable EE allows efficient evaluation of the system performance without the 
need of heavy Monte Carlo simulations. The simulation results demonstrate that the performance 
of the proposed scheme is almost identical to the near-optimal EE even for small system 
dimensions, with much reduced computational complexity. 

The rest of this paper is organized as follows: Section |n] describes a system model and the 
problem formulation. In Section Hill we present a low complexity beamforming design based on 
large system analysis, and simulation results are presented in Section |Vl Einally, in Section |Vll 
this paper is terminated with conclusions. 
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Fig. 1. The system model of M-cell MISO-IFBC 


Throughout the paper, we adopt upperease boldface letters for matrices and lowercase bold¬ 
face for vectors. The superscripts (■)^ and (•)'^ stand for transpose and conjugate transpose, 
respectively. In addition, || • ||, tr(-), [-j^ and [■]ij represent 2-norm, trace, the k-th element of a 
vector and the (i, j)-th entry of a matrix, respectively. Also, 1^ denotes an identity matrix of size 
d and 0^ means a zero matrix of size d x d. A set of N dimensional complex column vectors 
is defined by and |iS| indicates the cardinality of the set S. 

II. System Model 

In this paper, we consider an M-cell MISO-IFBC with bandwidth Wc where each BS equipped 

with Nt transmit antennas serves K users with a single antenna as shown in Fig. [T] Here user 

1 

{j,k) indicates the k-th user in cell j. Denoting hjnj,k = ^mZm,j,k where G is a 

deterministic transmit covariance matrix at the m-th BS and Zm,j,k ^ as the flat-fading 

channel vector from the m-th BS to the k-th user in cell j with the coherence time Tc, the 
received signal yj^k at user (j, k) is expressed by 

Vjjk ^j,j,k^j,k^j,k T ^ ^ ^rn,j,k^rn,n^m,n f^j^k 
{m,n)^{j,k) 



















5 


where Vm,n G equals the beamforming vector for user (m, n), Sm,n ~ CAA(0,1) stands for 

the complex data symbol intended for user (m, n), and rij^k ~ CA/”(0, represents the additive 
white Gaussian noise at user (j, k). Throughout the paper, we assume that the entries of ^m,j,k 
are uncorrelated Rayleigh fading according to CAA(0, em,j,k) where em,j,k indicates the pathloss 
from BS m to user (j, k). Note that our result can be easily extended to a more general channel 
model. It is also assumed that J2k=i order to satisfy per-BS power constraint 

Pr 

We assume single user detection at the receiver so that each receiver treats interference as the 
Gaussian noise. Thus, the individual rate of user (j, k) for given transmit beamforming vectors 
of all BSs {vm,n} is computed as 


log2 (1 T sinTj 


( 1 ) 


where sinrj^fc({vm,n}) represents the individual signal-to-interference-plus-noise ratio (SINK) for 
user (j, k) expressed as 


sinrj-fc({v,^,„}) = 



Kj,k^j,k\ 

2 

'n,{m,n)^{j,k) 

Vi H ■\r 



( 2 ) 


Then, the total amount of information transmitted during a time-frequency chunk TcWc is given 
by 


TcWc ^ log 2 (1 + si n rj-fc). [bits] 

Meanwhile, for designing an energy efficient transmission algorithm, we consider the power 
consumption model for BSs in [0. Thus, the total energy consumption during the time-frequency 
block TcWc is modeled as 

Pt = TcWc II ir + MNtPc + MPo j [Joule] (3) 

where ^ > 1 is a constant associated with the power amplifier inefficiency. Pc is defined by 
^ with Pc' being the constant circuit power consumption proportional to the number 
of radio frequency chains, and Pq denotes Pq = ^ with Pq accounting for the static power 
consumed at the BS which is independent of the number of transmit antennas. For example. 
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Pc includes power dissipation in the transmit filter, mixer, frequeney synthesizer, and digital-to- 
analog eonverter |l6l. 

Aeeordingly, the EE in bits/Joule is defined as the weighted sum-rate (WSR) divided by the 
amount of the energy eonsumption as 

/2({vj,fc}) C J2j,k + MNfPc + MPq 

Here, it is assumed that the positive weight term Wj^k is predetermined by a seheduler aeeording 
to the priority of user (j, k) [[5l. Thus, the EE maximization problem ean be formulated as 


/i {vyfc} 

max - - r- 

{v.'.fcl J2 


K 


s.t. ^ ||vj'fc||^ < Pj for j = 1, 


M. 


k=l 


(5) 


We notiee that problem ([5]) is in general non-eonvex due to eoupled interferenee among users 
and its fraetional form. Therefore, identifying a solution of this problem is quite eomplieated. 
As an alternative, by applying a transformation from the fraetional programming into LP and 
the WMMSE approaeh sequentially, a loeal optimal solution ean be obtained as in ll5l. 


III. Conventional Approach for Energy Efficiency 

In this seetion, we briefly review a eonventional approaeh for the EE maximization in Q, 
whieh employs a two-layer optimization strategy. Eirst, in the outer layer, the fraetional program¬ 
ming problem is transformed into a linear programming problem with a new parameter. Then, 
for the given parameter, the inner problem is solved by using the WMMSE method developed 
in ifT^ . Eventually, a final solution is found by inner and outer loops iterations. 

The optimization problem dS]) ean be transformed into a linear programming problem by 
introdueing a new parameter rj. Erom the relationship between the fraetional programming and 
the parametrie programming Q, the original problem ([5]) ean be reeast as the following equivalent 





max 


V 


( 6 ) 


{Vj.fe}, r^eM’*' 

S-t- /l({Vi,fc}) - vf2{{^j,k}) > 0 

K 

\\^j,kf < Pj for j = 1,..., M. 

k=l 

For a fixed value of r], we have a feasibility problem in {vj^} whieh eheeks if F{r]) > 0, 
where F{r]) indieates the optimal value of the following problem 

max -Vf2{{^j,k}) (V) 

K 

s.t. Y \\^j,k\\^ < Pj for i = 1 , ■ ■ ■, M. 

k=l 

From Theorem 1 in [|5||, F{r]) is shown to be a monotonically decreasing function with respect 
to T] and the equation F{r]) =0 has a unique solution. As a result, the optimal value of r] can 
be identified using one dimensional search algorithms such as a simple bisection method [fTOll . 

Next, the optimal beamforming needs to be determined in the inner problem for a fixed 77 . 
For the given 77 , the inner problem ([7]), excluding terms irrelevant to the optimization variables 
{vj,fc}, is rephrased as 


max G({vj-fc}) (8) 

K 

s-t- 

k=l 

where ^({vj-fc}) = T.j,kiRj,k - vChjA?)- 

Note that this problem is quite similar to the SR maximization problem except for the power 
term in the objective function. Thus, using the relationship between the SR and the WMMSE, 
a solution of the problem ([ 8 ]) can be computed from the following equivalent problem 

min Y Wj,k{ej,kSj,k - loga Sj.fc - 1) + vC\Wj,kf (9) 


where the mean square error Cj^k is given by 




Y1 

{m,n)^{j,k) 


|2 I |2 2 

1 + \Uj,k\ o- , 


( 10 ) 


and {sj,fc} and {uj^k} are auxiliary variables. 

The above problem is still non-convex in terms of {'Vj^k}, {sj,k}, {uj^k} jointly, making the 
direet optimization of the problem diffieult. However, sinee the problem is eonvex with respeet 
to eaeh of the optimization variables and {wyfc}, we ean solve the problem with 

one parameter by fixing the other two, i.e., the problem ean be calculated by alternating the 
optimization method. For given {vj fc}, the optimal Uj^k of the problem dH) is obtained by 


opt,'n _ 


^^m,n I I 


+ (T^ 


( 11 ) 


Furthermore, for fixed {wyfc} and {vj fc}, the optimal Sj ^ is expressed by 


. -1 


( 12 ) 


Then, once the values of {sj^k} and {wyfc} are given, the optimization of {vj fc} is decoupled 
among the BSs by substituting ^ in (fT^ into the objective function of the problem dH), and 
this leads to the following distributed optimization problems for the j-th BS 


K 


min Wj^kSj,k\Uj,kh.'^j^k^j,k - Ip + y^Wrn,nSrn,n\Urn,n\^\il^,rn,n^j,k\‘^ + VC\\^j,k\\‘^ |(13) 


k=l 

K 




k=l 


Denoting > 0 as the Lagrange multiplier corresponding to the power constraint, the first 
order optimality condition of the Lagrange function with respect to each ^ yields 

"^j,k ^j,k^j,kUj^k I ^ ^ '^m,n^m,n\Um,n\ ^j,m,rJ^j^rn,n i.VC T F'i)!- 1 (14) 

\{m,n) j 

where fij is chosen such that the complementary slackness condition of power constraint is 
fulfilled. Let Vj^kil^j) be the right-hand side of (fT4l) . If Yl!k=i — Pp *^ben = 0. 

Otherwise, can be found by using the bisection method which satisfies Ylik=i = 
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Pj. Therefore, a solution of for a given r] ean be eomputed by updating {s^^} 

and {wyfc} in an alternating fashion. 

In summary, a loeal optimal point of the EE ean be determined by two-layer optimization. 
However, the algorithm should be earried out in an iterative manner per each instantaneous 
channel realization by sharing global channel information among BSs. This leads to high com¬ 
putational complexity and signaling overhead. In the following, we will propose a new algorithm 
with low complexity and overhead which is more desirable in practical systems. 


IV. Proposed EE Scheme based on Earge System Analysis 


In this section, we propose an energy efficient scheme with low complexity in a finite 
dimension. After introducing conventional approaches, we describe the proposed method based 
on the asymptotic results of random matrix theory. Note that we consider the asymptotic regime 
where Nt ^ oo with ^ held at a fixed ratio to quantify beamforming parameters. The key 
idea is to combine large system analysis techniques with the WMMSE approach. Specifically, 
by applying the equivalence property between SR and WMMSE, the structure of the optimal 
beamforming vector is characterized with parameters related to beam direction and power. Then, 
by employing the asymptotic results of random matrix theory, the value of the parameters 
becomes deterministic which depends only on the second order channel statistics, and this leads 
to a significant reduction in the computational complexity compared to the system which utilizes 
instantaneous CSI. As will be shown later, the beamforming vector can be computed using only 
the optimized parameters and local CSI. 

Before explaining the algorithm, we provide useful results for solving the problem. Eirst, we 
identify the structure of the optimal beamforming based on in (fT4l) . Note that from 
the structure of beamforming can be paramterized with the power term pj ^ and the parameters 
related to the beam direction and \j as 


opt,ri 



(15) 
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where j3j^k and \j are given as 


m,n 

A, 


vC + f^j, 


(16) 

(17) 


and pj^k and Cj^k are denoted as = ||v°J’'||2 and cj^k = ^/\\{Y.{m,n) ^rn,niij,m,nh.'j^m,n + 
Ajl)“^hjfc||, respectively. Based on this structure, the normalized beam direction vectors ^ 
is defined as 

= Cj^k I ^ ^ (^^) 

\(m,n) / 

where l3m,n and Xj represent the parameters which control the leakage interference power level 
to other users adaptively. 

In order to quantify the component-wise impact on the performance for the given beamforming 
vectors we introduce the normalized channel gain matrix as 





lt.H |2 

0 ^1,1,2^1,1 

■ |hi,M,xVi,i|2 



|hl,l,lVi,2p 

ItiH v |2 

■ |hi,M,xVl,2p 

^ ^MKxMK 

(19) 


lt.H .Tr |2 





As mentioned before, in order to compute these instantaneous channel gains in a finite 
dimension, we exploit the results of the RMT for an asymptotic region. The ((j — 1) x J-f + 
k, {m — 1) X K + n)-th off-diagonal element of Q accounts for interference power at user (m, n) 
generated by the j-th BS for serving its supporting user (j, k), and the diagonal elements stand 
for the desired signal power. Employing the asymptotic results of random matrix theory, we 
arrive at the following lemma. 

Lemma 1: For fixed ratio with iV* —)■ oo, the deterministic equivalent of Q is obtained as 


Omk, 


g-g 


a.s. 


( 20 ) 
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where 


g° = 


Here, Dl^ and „ are given by 


n° 

^1,1 

T° 

-'1,1,1,2 

T° 



T° 

-'1,2,1,! 

-'^1,2 

T° 

^l,2,M,K 

^ ^MKxMK 

(21) 

T° 

'M,ir,i,i 

T° 

^M,K,1,2 

^M,K J 




D°u = 


[m 


o 'i 2 

'j,k) 


Nt^k ’ 


1°, 




j,k,m^n 


‘-j,k,m,n (IJ^R rn° 'l ’ 

Pm,nill'j^k,m,n) ^j,k 

where m°^ ^ ^ and ^lk,m,n are provided in the proof. 

Proof: See Appendix \M 

With g°, the deterministie equivalent of sinr^ ^ is derived as 


■ O Nt^cx) 

s\nrj^k - sinr ■ ^-)■ 0, 


where 




9j^k,j,kPj,k 


( 22 ) 

(23) 


(24) 


(25) 


^^{m,n)^{j,k) 9m,n,j,kP^ti^ ^ 

Here, 9^,1%,j,k represents the (m xn,jx k)-th element of g°. By the eontinuous mapping theorem 
m, one can show that the deterministic equivalents of the SR and the EE are expressed as 

9 j,k,j,kPjtk 


-^S '^j,k log2 [ ^ gO I 2 

Rl 


T] = 


CEj,kPj,k + MNtPc + MPo’ 


(26) 

(27) 


In what follows, based on the above asymptotic results, we optimize the EE performance in 
the large system limit instead of the original problem Q. The EE problem in the asymptotic 
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regime is formulated as 


max 7]° (28) 

K 

S-t- '^Pj,k<Pj Vj. 

k=l 


In the outer layer optimization, our algorithm is analogous to the eonventional approaeh shown 
in Seetion III. However, unlike the eonventional approaeh adopting a short-term strategy, we 
eonsider a long-term strategy in order to aehieve low eomplexity. Notiee that for the finite 
dimensional ease, the outer layer algorithm is required only to generate p and pass to the inner 
problem. Then, rj is updated aeeording to the feasibility, i.e., Fij]) > 0 or not, based on a 
solution of the inner problem. Similar to the finite ease, the feasibility in the asymptotie regime 
ean be determined whether F°{j]) > 0 or not, where F°{j]) is the optimal value of the following 
problem 


max 


s.t. 


Be - (c E P,.t + ^(iV,Po + a)) 

V j,k j J 

K 

Vj. 


k=l 


(29) 


Now, the only remaining work in the outer layer optimization is to identify the maximum 
value of r] for a biseetion method. It is elear that the maximum value ean be obtained as 

Eyfc'W^i,felog2(l + Sllhyi.fcP) 


Vmax 


EANtFc + Fo 


(30) 


Thus, applying Theorem 3.4 in [[T^ . it follows 


Vmax 


ZANtF, + Fo 


(31) 


where we have used — Nt 0. 

Next, we derive the optimal beamforming in the inner problem for a fixed p. Similar to ([8]), 
the objeetive funetion in the problem (l29l) ean be redueed to 


{Aj}, {pj,k}) = F.l- r/C 

j,k 


(32) 
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Then, by utilizing the WMMSE approach, the problem (1291) can be recast as 


/ 1 5^K',fc(ei,fc-log2Sj-fc-l) + r7CPi,fc) 


(33) 


j,k 

K 


S.t. 


J2P3^k<Pj Vj 


k=l 


where the mean square error ^ is given by 

9j^k,j,kPj,k ~ 

{m,n)^(j,k) 


2 o I 2 2 

^j,k9m,n,j,kP'rn,,n “T Uj j.(T . 


Thus, for any given {pj,k}, {/Sj.fc}, {Aj}, the optimal receiver filters of the problem (l33l) are 
obtained by 

y/ 9j^k,j,kPm,k 


yOPt - 
^3,k - 


9m,n,j,kP'kn,n T ^ 

Furthermore, the optimal Sj^k is expressed by 

1 


(34) 


„opi _ 

^ A U - 


(35) 


^-Uj,k^glk,j,kPj,k 

It is obvious that the optimal in (1331) is equal to from (fT^ for the finite case. 


Thus, we can calculate new Q° based on the updated }. 

Next, for given {syfc} and {uj^k}, the distributed problem for the j-th BS in the large system 
regime becomes 


{P: 


min ^ Wj,kSj,k\Uj,kJpj,k9lk,j,k “ iT + Wm,nSm,n\Um,n\'^g°m,n,j,k + VCPj,k (36) 

^ ^=1 \ (m,n)^(yfc) / 

K 

S.t. Pj,k < Pj- 


k=l 


Then, the optimal transmit power {pj^k} is written by 

'^j,kSj,kUj^k\/ 9j,k,j,k 


„opi _ 

Pj,k 


(37) 


'^m,n '^rn,nSrn,nU^ n9j^k,rn,n + 

Here, from (fTTl) . is equal to + where is determined by the complementary 

slackness of power constraint. 
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Algorithm 1 Outer Layer 

1 : Initialize r]min = 0 and T]max = Vmax in ([IB 

2: while |'?7maa; ^mm| ^ ^ tlO 

o. _ 'nmin~\~V max 

:>. If — 2 

4: Obtain the optimal solutions {Aj}, {pyfc} by Algorithmic 

5 : Compute F°{r]) 

6 : if F°{r]) < 0 then 

Vmax V 

8: else 

9- Vmin V 

10 : end if 

11 : end while 


Let us denote Pj,k{V'j) as the right-hand side of (iJTl) . If Ylk=i Pi,fc(0) < Pj, thenp^J = ppfc(O). 
Otherwise, we must have 


K 

= Pj- (38) 

k=l 

Aeeording to the monotonie property of the funetion <p(/ij) with respeet to pj, equation (1381) 
ean be efficiently solved by a bisection method. 


A. Overall Algorithm and Complexity Analysis 

The Algorithm [T] and |C describe the overall procedure of the proposed scheme. Here, 5 
indicates a predefined threshold. It is worth noting that our proposed algorithm depends only on 
the second channel statistics, and not on instantaneous channel realizations. More specifically, 
the conventional method determines {/Syfe}, {Aj} and {pj^k} per each channel realization, while 
the proposed algorithm does so only when the second order statistics changes, i.e., signal- 
to-noise ratio (SNR) changes. Once the optimal and {pj,k} are determined, we 

can construct the beamforming vectors based only on local CSI without additional complexity. 
Thus, the proposed algorithm dramatically reduces the computational complexity compared to 
the conventional scheme. 

In what follows, we compare the complexity of the conventional scheme with that of the 
proposed method. For comparison, the overall computational complexity can be characterized by 
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Algorithm 2 Inner Layer 


1 : Initialize n = 1, = Wj^k, = vC,P^ik = ^ 

2 : SetG"({f)j5},{A<“'},{pf>}) = 0 
3 : Compute Q° 

4 : While |G“({/?‘"J}. {Af}. {p<”i}) - G“({/3<”-'>}. {A^”-*'}. {p^”-*'})! > i do 
n n + 1 


5: 

6 : 

7: 

8 : 

9: 


Update with (1^ and (1^ for all j, k 

Caleulate 

Compute with 


for all j, /c 


10 : end while 


the multiplication of the following terms: the execution rate of the algorithm, the iteration number 
of the outer layer, the iteration number of the inner layer, and the complexity of inner layer 
optimization per each iteration. Here, the per-iteration complexity of the outer layer optimization 
is ignored since the calculation using the bisection method is relatively simple. 

First, the per-iteration computational complexities of the conventional scheme and the proposed 
method are 0{M‘^K‘^Nf) and 0{M‘^K'^), respectively. The difference between the two schemes 
comes from an inverse operation of an Nt x Nt matrix. The conventional algorithm requires 
the inverse operation to generate the beamforming vectors {v^ in (fT^ . while the proposed 
scheme does not need as can be seen in Algorithm [21 Thus, a complexity gain becomes larger as 
Nt increases. Next, for the inner and outer layer algorithm, the required number of iterations of 
both schemes are quite similar in average sense, since both of them are based on the bisection 
method and the WMMSE algorithm. 

Compared to the conventional scheme, the factor which reduces the complexity the most in 
the proposed algorithm is the execution rate of the overall algorithm. The update rate of the 
proposed algorithm depends on how often the second order channel statistics changes, and thus 
is much slower than that of the conventional method which needs to update at each realization. 
This is because large scale fading varies with tens of seconds, while small scale fading changes 
with few milliseconds in general wireless environments [|20||. Thus, the coherence time of small 
scale fading is typically 1000 times smaller than that of large scale fading. As a result, for 
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TABLE I 
System setup 


System bandwidth (W) 

20 MHz 

The number of user drops 

10 

The number of channel realizations per user drop 

100 

The number of Tx antennas for each BS Nt 

4 

Cell radius R 

500 m 

Minimum distance from BS to each user Rmin 

35 m 

Pass loss exponent a 

3.8 

Transmit power constraint per BS Pj 

26 ~ 46 dBm 

Circuit power per antenna Pc 

30 dBm 

Basic power consumed at BS Pq 

40 dBm 

Noise figure (Np) 

7 dB 

Noise power 

-94 dBm 

Inefficiency of the power amplifier ( 

2 


example, with M = 3, Nt = A, and K = 3, the CPU running time of the conventional 
EE algorithm is about 300 times more than that of the proposed scheme. Therefore, we can 
verify that the proposed algorithm greatly reduces the computational complexity compared to 
the conventional scheme. We shall show in the simulation section that our proposed algorithm 
exhibits the performance almost identical to that of the conventional algorithm while requiring 
significantly reduced complexity. 


V. Numerical results 

In this section, we evaluate the EE performance of the proposed beamforming scheme. We 
consider a cooperative cluster of M = 3 hexagonal cells for Monte Carlo simulations. These 
simulations are carried out with the parameters listed in Table HI unless specified otherwise. The 
pathloss from BS m to user (j, k) em,j,k is given as 10 log^o ^m,j,k = —38 log^o dm,j,k — 34.5 in 
decibels, where dm,j,k in nieter indicates the distance from BS m to user (j, k). Also, the noise 
power can be calculated as = —174 + 10 log(>V) + Np in dBm where W means the system 
bandwidth and Np denotes the noise figure. 

Eirst, we illustrate the convergence of the proposed algorithm. Eor the case of the outer layer 
optimization, the optimal r] can be found based on one dimensional line search without loss of 
optimality, and thus the convergence is guaranteed. On the other hand, the inner layer algorithm 
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20 


10 15 
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Fig. 2. Convergence examples of the inner layer algorithm in the proposed scheme 


cannot achieve the global optimal value due to the non-eonvexity of the problem. However, the 
eonvergenee to a loeal optimal point is guaranteed by virtue of the WMMSE approaeh ETI . Fig. 
[2] plots the objeetive funetion G° in (l32l) with respeet to the number of iterations for K = A. The 
eonvergenee trend varies with parameters sueh as power, user position and 7]. In this figure, the 
eurves eorresponding to 5 different user drop events are plotted by fixing eertain r] and Pj = 46 
dBm. As shown in this plot, the inner layer algorithm eonverges to a stable point with about 10 
iterations. 

Fig. [3] exhibits the average FF performanee of various beamforming sehemes as a funetion 
of Pj for K = 3 with te = [1 2 3]. For eomparison, we first present the FF performanee of the 
following beamforming sehemes. 

• Maximal ratio transmission (MRT): the beamformers are aligned with the eorresponding 
ehannels. 

• Zero-foreing beamforming (ZFBF): the signal to unintended users is nullified. 

• Conventional virtual SINK (VSINR): the VSINR is maximized with non-weighted eoeffi- 
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Fig. 3. Average EE performance of various beamforming strategies with K = 3 


dents, i.e., (3j^k = 1 and Xj = ^ for all i,j l[22]l . Il2^ . 

. WMMSE algorithm: beamformers are designed to maximize the SR by using the WMMSE 
approaeh ifT^ . 

• Conventional EE algorithm: the algorithm based on the WMMSE approaeh is adopted to 
maximize EE as deseribed in Seetion III. 

. Proposed EE algorithm: the proposed algorithm performs with adaptive control of {Aj} 

and {pj^k} based on second order channel statistics for the EE maximization. 

Eor both the WMMSE and the conventional EE schemes which achieve the local optimal solution, 
a solution of the VSINR scheme is adopted as an initial point. Surprisingly, we can see that 
the proposed scheme achieves near-optimal performance with much reduced complexity for all 
simulated transmit power constraint ranges. It is emphasized again that our proposed algorithm 
is performed only when second order channel statistics changes and the constant values of 
fc}, {Aj} and {pj,k} are employed for generating beamforming vectors as long as the statistics 
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Fig. 4. The EE performance for different number of users 


remains unchanged. This results in a significant computational complexity reduction compared to 
the conventional EE scheme which should be carried out in every channel realizations. Moreover, 
we can observe the trade-off relationship between the performance and the complexity for 
various beamforming strategies. Simple beamforming schemes such as MRT, ZEBE, and VSINR 
require lower computational complexity to comprise the beamforming structure than the proposed 
scheme but with poor performance. In the ZEBE case, the EE performance is mainly degraded 
by the deficiency of dimension for nullifying the unintended user signals. We also observe that 
the WMMSE schemes designed for spectral efficiency maximization produce much worse EE 
performance compared to the proposed scheme. Especially, a performance gain of the proposed 
EE scheme is about 209% at Pj = 46 dBm Vj. 

Also, the average EE performance for the conventional and proposed EE algorithms is illus¬ 
trated for various number of users K in Eig. HI We can see that the EE performance of these 
two algorithms increases with the number of users. Moreover, the performance of the proposed 
EE algorithm has a small gap compared to the conventional algorithm, which is less than 4% 





















































20 



Fig. 5. The EE performance with transmit antenna correlation (K = 2) 

for all cases. 

In Fig. [51 we demonstrate the EE performanee of the proposed seheme for the eorrelated 
transmit antenna ease. The transmit eovarianee matrix is set as the exponential eorrelation model 
whieh is given by [R]ij = pi*”-’! with i, j = 1, • • • , iVt and p G [0,1). In this eorrelation model, 
the average EE performanee is enhaneed when the eorrelation eoeffieient p grows at a high 
transmit power region. On the eontrary, the opposite trend is observed at a low transmit power 
region. This is due to the faet that the EE performanee ean be mainly affeeted by an array gain 
at the low power region. Eor the high BS power region, the effeet of spatial multiplexing takes 
a key role of the EE performanee. 

In what follows, we validate the aeeuraey of the deterministie equivalent of EE eompared to 
true EE. Eig.j^eompares the average EE with the deterministie approximation for Nt = 10, 20 
and 40 with the fixed ratio ^ = 2. In this plot, eaeh eurve eorresponds to a partieular drop of 
users for eaeh Nt. The error bars indieate the standard deviation of the simulation results. As 
shown in this plot, the deterministie equivalent of EE provides a very aceurate approximation. 
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Fig. 6. Comparison between the average EE and the deterministic equivalent for -^ = 2 

It can be seen that the approximation lies within one standard deviation of the Monte Carlo 
simulations and the standard deviation becomes smaller as Nt increases. Also, we can check that 
the maximum value of EE gets larger as Nt grows. This comes from an increased multiplexing 
gain, (i.e. pre-log term) as Nt and K grow larger. Erom the plots, it is observed that the EE 
performance curve increases up to a certain point, and after that it is saturated. Investigation of 
the saturation point will be an interesting future work. 

VI. Conclusions 

In this paper, we have proposed a low complexity beamforming scheme for MISO-IEBC. 
With the parametrization of the beamforming vectors by the scalar values, we have found the 
optimal parameters to maximize the EE in the asymptotic regime. Our solutions depend only on 
second order channel statistics, not on instantaneous CSI, and thus the parameters are computed 
only when channel statistics changes. As a result, the computational complexity is significantly 
reduced compared to the conventional method. Through simulations, we have confirmed that the 
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proposed schemes with the asymptotic results provide the near-optimal EE performance even 
for the finite system case. Additionally, the proposed scheme allows efficient calculation of the 
system performance without resorting to heavy Monte Carlo simulations. 


Appendix A 
Proof of Eemma[I] 

We will derive the deterministic equivalents of the desired signal power and the 

interference power subsequently as in ifTTll . Eor simplicity, we assume Rj = R. 

1) Deterministic equivalent for Eor given {/Sj^k} and {A^}, is written 

by 






.k^jM 


(Aj + AjI) 


j,k\ 


|h,-.>(A, + A,I)"^h,,,,p 

(Aj + AjI)"^ 

K,k{^jk + \,ir^h,,^k\^ 


(Ajfc + AjI) 


(39) 


where A^ = ^ ^jk = J2{m,n)^(j,k) and the last equality 

comes from the Sherman-Morrison matrix inversion lemma. 

Eirst, applying Theorem 3.4 in lfT9l to the term {Ajk + Ajl)~^ ^j,j,k in the numerator of 
(l39l) yields 


(A,» + hjj,t - (REji) 


(40) 


where = 
follows 




0m,n 


hym,nh^^„ + ^l) . By employing Theorem 1 in (SI, it 


Nt Nt' 


-1 


^j,j,k (Ajfc + AjI) i^j,j,k 


d.s. ^ 
m, j, —)■ 0 




(41) 


where fctr (^R0(£jfc, ^) j and Cjk = -Cj\{ejj^kl3j,k} with Cj = {e^-e^-i, 2 /Si, 2 , • • • 

j^j,M,Kl3M,K}- Here, (j){S,p) is defined as 


\ Si^S 


1 + Ci 


-1 


pi 


(42) 








23 


where S equals a set with non-negative elements s* for i = 1,..., |iS|, p represents a positive 
sealar value and e/s are unique positive solutions of the fixed-point equations 


Ci = Sitr (R0(5,p)). 


Next, for the denominator in (1391) . Theorem 1 in [|24ll leads to 

Ki.k (Ait + h,,,,s - ^tr (RSJJ ^ 0. 


Then, adopting Theorem 2 in If2^ . we ean write 


h,- (A,, + A,I)-^ h,,,, - ^ 0 

where '^°jk = {^jk, ^))- Here, (j)'{S, p) is denoted as 


'>'(S,p) = <j>(S,p) I 1+ -A I <l>(s,p) 


N. (1 + 


and e' = [e [,..., is expressed by 


e' = (l|5|-J) V, 

where J and v are eomputed as 

_ 5igjtr(R0(5,p)R0(>S,p)) 

^ Nt{l + ejy 


for i,j = 1,..., |5| 


= [sitr (R0(5,p)^) ,...,S|5|tr(R0(5,p)^)]’^. 


Finally, eombining (|431) and (|4TI) . we have 


IhH 12 




(mlkf a.s. 


Nt J,k 


(43) 


(44) 


(45) 


(46) 


(47) 
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2) Deterministic equivalent for The interference term can be written 


as 


IV.H V |2 - 
j,k\ 


|h"„^„(Aj + A,I)-‘hj.j.tP 
II (A, + A^I)-‘ 

^j,m,n {-^jk + -^jl) ^j,j,k^jj,k {-^jk + -^jl) i^j,m,n 
^j,j,k i-^jk + -^jl) ^j,j,k 

^j,m,n {-^jkmn + -^jl) {-^jkmn + -^jl) ^j,m,n 


(48) 


(l + l^m.,vih-j^rn,n {-^jkmn + -^jl) (-^jk + ^j,j,k) 

where Ajkmn = T.{i,q)^{j,k),{m,n) second and third equality come from the 

Sherman-Morrison matrix inversion lemma with respect to h^j ^ and hj^rn,n, respectively. 

First, by applying Theorem 3.4 in lfT9]l to the numerator in (l48l) twice with respect to hj j ^ 
and hym,n> we can easily show that 

h"™,„ + A,I)-‘ (A,t„„ + A,I)-' 0 (49) 


where 12^kran is defined as 'Sjkmn=\^jf^Nt+Yl 
in lf24]l . it follows 


.i^,q)^(j,k),(m,n) ^K^,q^l^,q) • From Thcorcm 1 


-1 


N? 


tr - -»» 


k.m.n 


(50) 


where ^j^k,m,n ^iqtk^j,m,nX'^ ^ ^^jkmnt JVf)^ with C^jkmn C,jk\\^ej^rn,,ril^m,,n\- 

Similarly, in the denominator of (l48l) . one can show that 


Fj,m,n {Ajkmn + AjI) hym,n ^j,k, 


(51) 


where ^ 


converges almost surely to Then, combining all results generates 


IViH Va. |2_ 
j,k\ 


)). Also, from (1451) . we know that the second term 


(52) 


4/°i. 

j,k^m^n 




^0, 


and this concludes the proof. 
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